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ABSTRACT 

We present results of radial-velocity follow-up observations for the two Kepler evolved stars Kepler- 
91 (KOI-2133) and KOI-1894, which had been announced as candidates to host transiting giant plan¬ 
ets, with the Subaru 8.2ni telescope and the High Dispersion Spectrograph (HDS). By global modeling 
of the high-precision radial-velocity data taken with Subaru/HDS and photometric ones taken by Ke¬ 
pler mission taking account of orbital brightness modulations (ellipsoidal variations, reflected/emitted 
light, etc.) of the host stars, we independently confirmed that Kepler-91 hosts a transiting planet 
with a mass of 0.66Mjup (Kepler-91b), and newly detected an offset of ~20 m s“^ between the radial 
velocities taken at ^ 1-yr interval, suggesting the existence of additional companion in the system. 

As for KOI-1894, we detected possible phased variations in the radial velocities and light curves with 
2-3(7 confidence level which could be explained as a reflex motion and ellipsoidal variation of the star 
caused by the transiting sub-saturn-mass 0.18Mjup) planet. 

Subject headings: stars: individual: Kepler-91 (KOI-2133) — stars: individual: KOI-1894 — planetary 
systems — techniques: radial velocities 


1. INTRODUCTION 

Detecting planetary transits is highly valuable not 
only because it makes an independent confirmation of 
a planet from radial-velocity observations but also be¬ 
cause the photometric transits provide unambiguous in¬ 
formation on planet masses and radii, and th us mean 
density and interior structure (e.g. iWinn I I 2 OIIL and ref¬ 
erences therein). Various high-precision follow-up studies 
for transiting planets can also uncover planetary atmo¬ 
spheres and the (mis-)alignment between the stellar spin 
and planet ary orbital axes via the Rossiter-McL aughlin 
effect (e.g. iSeaeer fc Deming|l201^ IWinn 1 1201 111 . These 
properties of transiting planets provide valuable hints for 
planet formation and evolution processes including dy¬ 
namical interaction between planets. 

Although transiting planets are as such important, 
it is difficult to detect them around evolved stars, es¬ 
pecially giant stars, due to the large sizes of the host 
stars. Targets for on-going radial-velocity surveys for 
planets around giaii t s have typical radii ( jf ^10 Rrrj (e.g. , 
da Silva et al.ll20(Ifil iTakeda et al.l 120081 : iLiu et all 1201(11 : 
Wang et al.ll2011l : lZielinski et al.ll2012f) . and then the rel¬ 
ative flux variation of such a giant host star caused by a 
transit of a Jupiter-sized planet is only ^ 1 x 10“"^, which 
is comparable to that for transit of an Earth-sized planet 
across a solar-type star. It is impossible to detect such a 
transit from ground-based photometry and thus no tran¬ 
siting planets had been found around giant stars. Our 
understanding of properties of planets around such stars 
are thus far behind from those around solar-like stars, 
although tens of planets have been found around giants 
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by precise radial-velocity surveys (e.g. iSato et aI1l20I3l 
and references therein). 

The Kepler mission was successful in detecting tran¬ 
siting planets with very high photometric precision (^ 

2 X 10“®) since 200!fl. To date it has discovered more 
than 4000 transiting planet candidates from sub-Earth- 
size to super-Jupit er-size in ^0.3-2000 d orbits (e.g. 
iBatalha et al.l I20I3I1 . Although most of them are or¬ 
biting solar-type stars, Kepler^s photometric precision is 
high enough to detect a Jupiter-sized planet transiting a 
giant star. 

As expected, Kepler has identified several planet can- 
didates around gian t stars with radii larger than 5 Rq 
(jBatalha et al.l 12(11^ . The planet candidates have radii 
comparable to or larger than Jupiter’s, and it is in¬ 
teresting that many of them have short-period orbits. 

Since such short-period planets have rarely been found 
by radial-velc)city s urveys ar(3und evolved stars (e.g . 
[Johnson et al.|[200^ ISato et ^120081: 1 Jones et'^l2014h . 

they will provide us unique opportunities to investigate 
planet formation and evolution processes around such 
stars. 

Here we report the results of radial-velocity follow¬ 
up observations using Subaru 8.2m telescope for two 
of the Kepler evolved stars, Kepler-91 (KOI-2133) and 
KOI-1894. Thanks to the high precision in their radial- 
velocity measurements, we independently confirmed 
a jovian planet (Kepler-91b, KOI-2133.01) previous ly 
reported around K epler-91 (|Lillo-Box et al.l l2014al la: 
[Barclay et al.l[20T3l and newly found a hint for the ex¬ 
istence of additional companion in the system. We also 
detected a possible sub-saturn-mass planet around KOI- 
1894 (KOI-1894.01) with 2-3 (t conficience level. 

The rest of the paper is organized as follows. The 
adopted stellar parameters for the two stars are presented 
in section [21 The observations are described in section [3] 

^ NASA Exoplanet Archive (http://exoplanetarchive.ipac.caltech.edu/ ) 
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and the results of global analysis for the radial velocities 
and light curves are presented in section |4l Section [5] is 
devoted to discussion and summary. 

2 . TARGETS 

Kepler-91 (KIC 8219268, KOl-2133; Kp=l2.A9®) 
was identified as a candidate star hosting a tran¬ 
siting short-period jupiter-sized planet (KOl-2133.01; 
Rr, = 18.24 Rm P = 6.2465798 ± 0.0000821 d) 
by iBatalha et all (I2013D . After that, iLillo-Box et al'l 
(j2014ar) reported the confirmation of its planetary na¬ 
ture based on the detailed analysis of orbital brightness 
modulation seen in the light curve caused by ellipsoidal 
variation, Doppler boosting, and ref l ected / emitted light 
from planet (e.g. iFaigler fc Mazehl 120111 : iMazeh et al.l 
1201211 . and then the planet was named Kepler-9 lb 
with the radius Rp = 1.384lg;Q54 Rjup and mass 
Mp = 0 . 88 lo ;33 Mjup. Although lEsteves et al.l (l2013f l 
and iSliski fc Kipping! (1201411 cl aimed a possible non - 
plan ertary nature for the system. iLillo-Box et al.l (j2014bf l 
and [Barclay et al.l (|2014f l very recently reconfirmed the 
planetary nature via radial-velocity measurements with 
a precision of ^100 m s“^ and ^20 m s“^, respectively, 
and obtained the planetary mass of 1.09±0.20Mjup and 
0.73±0.13Mjup, respectively. 

The stellar parameters (effective temperature Teff, sur¬ 
face gravity log q, mass Mj,, and radiu s i?*) of Kepler-91 
were reported bv IBatalha et al.l (|2013[1 to be = 4712 
K, logg = 2.85 cgs, M* = 2.25 Mq, and i?* = 9.30 Rq. 
After that, the parameters have been up dated by the 
spectro scopic and asteroseismic anal yses in iHuber et al.l 
(l2013all and iLillo-Box et al.l (|2014all . which are consis¬ 
tent wit h each other. Thus, we here adopted the values 
listed in iLillo-Box et al.l (j2014all : = 4550 ± 75 K, 

logt; = 2.953 ± 0.007 cgs, M* = 1.31 ± 0.10 Mq, and 
i?* = 6.30 ±0.16 Rq. 

KOl-1894 (KIC 11673802; 74:^=13.427) was also re¬ 
ported to be a planet-host candidate having a tran¬ 
siting short-period jupiter-sized planet (KOl-1894.01; 
Rn = 16.29 fim, P = 5.2880157 ± 0.0000447 d) by 
IBatalha et al.l (j2013ll . though radial-velocity follow-up 
observations for the star have not been reported yet. 
The stellar p a ramet ers of KOT1894 were derived by 
IBatalha et al.l (1201311 to be Teg = 4815 K, log 5 = 2.87 
cgs, M* = 2.02 Mq, and i?* = 8.62 Rq, and they 
have been upd ated by the spectrosc opic and asteroseis¬ 
mic analyses in iHnber et al.l (l2013a[l as Teg = 4992 ± 75 
K, M* = 1.410±0.214 Mq, and A* = 3.790±0.190 Rq. 
Here we adopted the up dated values in this paper. 
Recently iLaw et all (1201411 reported a non-detection of 
blended stars for KOT1894, which could have been phys¬ 
ically associated companions and/or responsible for tran¬ 
sit false positives if they were within ^0.15"-2.5" sepa¬ 
ration and with magnitude difference up to Am ~ 6 , by 
high-angular-resolution AO imaging. The parameters for 
the two stars we adopt here are summarized in Table [T] 

3. PRECISE RADIAL VELOCITY 
MEASUREMENTS WITH SUBARU/HDS 

We obtained high-precision radial-velocity data for the 
stars with the 8.2m Subaru telescope and the High Dis- 

® Kepler mag 


TABLE 1 

Stellar Parameters for Kepler-91 and 
KOl-1894 


Parameter 

Kepler-91 

KOl-1894 

TefI (K) 

4550±75^ 

4992±75‘> 

logs (cgs) 

2.953±0.007=>- 

2.87'= 

M* {Mq) 

1.31±0.10=‘ 

1.410±0.214‘' 

R* (Rq) 

6.30±0.16^ 

3.790±0.190‘' 


^ ILillo-Box et al.l J2014al l 
^JHuber^t_ayj 2013^ 

^ fBataIha*et~S (I2013l ) 


persion Spectrograph ("HDSi lNoguchi et'mi2002[ l in 2013 
and 2014. We used the setups of Stdl2b (2013 June 28-30 
and 2014 July 8, 10, 14-16) and Stdl2a (2013 December 6 
and 2014 July 12-13), which simultaneously cover a wave¬ 
length region of 3500-6200A and 4900-7600A respec- 
tively, the image sheer No.2 (134^2: iTaiitsu et 'n]l2012ll 
yielding a spectral resolution {R = A/AA) of 80000, and 
an iodine absorption cell (I 2 cell; iKambe et al.l 1200211 
for precise radial-velocity measurements. We obtained 
a total of 29 and 18 data points for Kepler-91 and 
KOl-1894 with typical signal-to-noise ratio of S/N=30- 
80 pix“^ and 30-45 pix“^, respectively, by an exposure 
time of 1200 sec depending on weather condition. The 
reduction of echelle data (i.e. bias subtraction, flat- 
fielding, scattered-light subtraction, and spectrum ex¬ 
traction) was performed using the IRAJ0 software pack¬ 
age in the standard manner. 

We performed radial-velocity analysis for I 2 - 
superposed stellar spec t ra (s tarTp ) by the metho d 
described in iSato et all (1200211 and iSato et al.l (1201211 . 
which is based on th e method bv iButler et al.l (jl996ll and 
iValenti et al.l (119951 1. A star±l 2 spectrum is modeled as 
a product of a high resolution I 2 and a stellar template 
spectrum convolved with a modeled instrumental profile 
(IP) of the spectrograph. We obtain the stellar template 
spectrum by deconvolving a pure stellar spectrum with 
the IP estimated from an U-superposed Flat spectrum. 
We achieved a radial-velocity precision of 5-17 m s“^ 
for Kepler-91 and 9-15 m s“^ for KOl-1894. The 
derived radial velocities are listed in Table [5] and Table 
[3] together with the estimated uncertainties, and are 
plotted in Figure [T] and Figure HI 

4. RADIAL VELOCITY AND LIGHT CURVE 
ANALYSIS 

4.1. Method and Light Curve Reduction 

The variations in the observed radial velocities for both 
targets show a sign of planetary companions to Kepler- 
91 and KOl-1894. In order to obtain accurate and pre¬ 
cise estimates for system parameters of those systems, we 
here present a global analysis that makes use of all the 
available information from the Kepler photometry and 
our spectroscopy. 

As is well known, a very precise light curve of a star 
orbited by planet (s) shows a periodic modulation due 

® IRAF is distributed by the National Optical Astronomy Ob¬ 
servatories, which is operated by the Association of Universities 
for Research in Astronomy, Inc. under cooperative agreement with 
the National Science Foundation, USA. 
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TABLE 2 

Radial Velocities for Kepler-91 


BJD-2450000 

Velocity (m s ^) 

Uncertainty (ms 

6473.00247 

43.20 

4.61 

6473.01719 

33.85 

4.58 

6473.03192 

45.26 

5.03 

6473.96075 

17.35 

7.36 

6473.97549 

8.93 

7.34 

6473.99022 

30.04 

6.39 

6474.95764 

-24.28 

5.43 

6474.97235 

-22.62 

5.03 

6474.98708 

-39.46 

5.50 

6633.69753 

-42.12 

17.38 

6847.93547 

42.30 

10.56 

6847.95020 

37.46 

13.63 

6847.96513 

11.49 

8.76 

6847.98006 

29.77 

7.96 

6847.99479 

17.07 

6.12 

6848.00952 

19.94 

8.61 

6848.08121 

23.25 

6.59 

6849.75412 

-52.12 

5.23 

6849.76886 

-36.10 

5.06 

6849.78360 

-34.47 

5.48 

6853.04854 

-40.23 

13.48 

6853.97273 

5.73 

5.72 

6853.98746 

9.02 

4.73 

6854.00219 

7.66 

5.17 

6854.84234 

2.14 

4.99 

6854.85708 

0.36 

5.25 

6854.87181 

-6.69 

6.03 

6855.80523 

-43.37 

6.63 

6855.81828 

-41.91 

5.28 


TABLE 3 

Radial Velocities for KOI-1894 

BJD—2450000 Velocity (ms“^) Uncertainty (ms“^) 

6473.04710 

-9.95 

9.45 

6473.06182 

-0.82 

11.01 

6473.07656 

-11.81 

9.94 

6474.00523 

-5.31 

10.97 

6474.01995 

1.95 

10.25 

6474.03468 

17.33 

10.36 

6474.04941 

5.73 

10.49 

6475.00210 

-3.24 

10.03 

6475.01683 

-21.07 

9.68 

6475.03159 

-13.24 

8.79 

6475.04632 

-6.92 

9.62 

6848.09625 

-4.14 

15.02 

6849.94269 

0.52 

10.34 

6849.95743 

13.69 

9.86 

6849.97216 

10.09 

9.04 

6855.88939 

0.65 

9.63 

6855.90411 

-3.06 

8.95 

6855.91885 

6.48 

9.09 


to several astrophysical effects: the ellipsoidal variation, 
Doppler boo sting, and reflected/emission light from th e 
planet le.g.. iFaigler fc Mazehll201ll : iMazeh et alJ[^ni2[l . 
The last two effects (Doppler boosting and planetary re¬ 
flection/emission) are synchronous with the planet’s or¬ 
bital period Porb, with different peak locations along the 
planet phase {(j) ^ 0.25 for boosting and (j) ^ 0.5 for 
the planetary reflection/emission, respectively). On the 
other hand, the ellipsoidal variation, which is caused by 
a tidal distortion by the planet’s gravity, have two flux 
peaks {(j) ^ 0.25 and 0.75) within one orbit. Thanks to 
the different phase dependence of these three effects, a 
very precise light curve and its modeling enable us to dis¬ 


tinguish these three effects, and we can extract physical 
parameters such as the planet-to-star mass ratio q and 
scaled semi-major axis a/R^,, where a is semi-major axis 
and i?* is stellar radius. Moreover, for transiting systems 
as in the present cases, incorporating the transit and/or 
secondary eclipse model into the above phase-curve vari¬ 
ation lets us learn more about the planet properties (e.g., 
planet-to-star radius ratio and orbital inclination 

G) ■ 

To obtain the phase-folded light curves for Kepler- 
91 and KOT1894, and estimate system parameters, 
we downloaded all available public light curves (QO- 
Q17) from Kepler MAST archive. While only long- 
cadence data were available for Kepler-91, KOT1894’s 
light curves involve some short-cadence data. Adopt¬ 
ing the PDC-SAP flux data, for which unphysical ar¬ 
tifacts are detrended, we reduced the light curves by 
the following procedure. First, after removing plane¬ 
tary transits, we further detrended and normalized the 
light curve for each quarter by fitting it with a fifth- 
order polynomial so as to remove the long-term trends 
that were not removed in the PDC-SAP flusQ. This 
process was repeated implementing a 5cr clipping to re¬ 
move outliers. We then combined the light curves for 
all quarters and phase-folded them with the ephemerides 
derived by the official Kepler team using the Q1-Q16 
data; orbital period Porb = 6.24668005 (±0.00002647) d 

and transit center = 2454969.38661 (±0.00346) BJD 
for Kepler-91, and Porb = 5.28789787 (±0.00001198) d 

and ri°^ = 2454968.36353 (±0.00189) BJD for KOI- 
1894 (from NASA Exoplanet Archive). We later consider 
the impact of incorrect ephemeris. Finally, the folded 
light curves were binned into 300 and 250 phase bins for 
Kepler-91 and KOI-1894, respectively. These bin num¬ 
bers were adopted so that each bin approximately covers 
a cycle span of long-cadence data (^ 30 minutes). The 
flux error for each bin ctlc j which is the standard devia¬ 
tion of the mean flux, was computed based on the disper¬ 
sion of the flux values within the bin. The long-cadence 
and short-cadence flux data were separately folded and 
binned for KOI-1894, but the binned light curve of the 
short-cadence data was much noisier than that for long- 
cadence. We thus decided to ignore the short-cadence 
data in the following analysis. 

The analysis be low is based on the method by 
iHirano et al.l (1201511 . who employed the EVIL-MC model 
(| Jackson et al.ll201^ for the phase-curve variation (i.e., 
ellipsoidal variations, Doppler boosting, and planetary 
reflection) with some revisions (e.g., application to an 
eccentric orbit). This phase-cu rve model i s mult iplied 
by the analytic transit model bv lOhta et al.l (|2009(1 . and 
the relevant parameters (e.g., a/P* and io) are simulta¬ 
neously determined. In addition, we also model and fit 
the observed radial velocities Uobs- Thus, the statistics 
in the present case is 

/,,(*) _ „(*) '12 ( Ai) _ Ai) ',2 

2 _ \ ^ V^obs ^model/ , \ ^ w obs J model' /-i \ 

^ ~ {i)2 ^ {i)2 ’ 

i ^RV i ^LC 

^ Note that since the orbital periods of our targets are much 
shorter than the time span of one quarter, the flux modulations due 
to ellipsoidal variations, Doppler boosting, and planetary reflection 
are retained in this process. 
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Fig. 1.— Phased radial velocities for Kepler-91 in the absence 
of radial-velocity trend 7 . The data in 2013 and 2014 are shown 
by green and red points, respectively. The solid line is the best-fit 
radial-velocity curve to the data obtained by global modeling of 
photometric and radial-velocity data without considering 7 . 


where and tr^y are the i-th observed radial velocity 

and its error, and and ct^q are the i-th observed 
light curve flux and its error, which are the mean flux 
and i ts standard deviati on in the f-th bin, respectively 
(e.g., iHirano et al.ll20lj1l . The radial velocity is modeled 
as 


Wmodel = .ftl{cOs(/ + w) + ecosw} + 7, (2) 

where K, /, e, w, and 7 are the radial velocity semi¬ 
amplitude, true anomaly, orbital eccentricity, argument 
of periastron, and radial velocity offset of our dataset, 
respectively. F or the light c u rve m odel /model, we refer 
the readers to IHirano et al.l (|2015ll . for details. In the 
next subsections, we describe the fitting procedure for 
each of the targets. 


4.2. Kepler-91 

The binned light curve of Kepler-91 in Figure[5]shows a 
clear pattern of phase-curve variation; the double peaks 
at (j) ^ 0.25 and 0.75 are representative of the ellipsoidal 
variation. In addition. Figure [H suggests a possible de¬ 
tection of secondary eclipse at ^ 0.5, which means the 
reflected/emitted light from the planet is visible in the 
folded light curve. The r efore, following iJackson et al.l 
(|2012fl and IHirano et al.l (|2015h . we model the planet’s 
light Fp by the following expression: 


1 -I- e cos / 


2 

sin(/ -I- oj 


2nA^), (3) 


1.0002 

1.0001 

1 

0.9999 

0.9998 

0.9997 

0.9996 

0.9995 

0.9994 

0.0001 

5e-05 

0 

-5e-05 

- 0.0001 
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Fig. 2. — Phased light curve for Kepler-91 (red points). The solid 
line is the best-fit light-curve model to the data obtained by global 
modeling of photometric and radial-velocity data considering the 
velocity trend (7 = —0.0612 m s“^day“^). 


where Fi is the flux variation amplitude of planet’s re¬ 
flected/emitted light and Fq is the planetary flux off¬ 
set ari sing from the homoge n eous surface emission . Fol- 
lowing iDemorv et al.l (120131) . lEsteves et al.l (|2014D . and 
iFaigler fc Mazehl (|2014l I), we here introduce the “phase- 
delay” Acj) of the brightest part on the planet surface 
from the substellar point. The planetary light Fp is 
added to the flux model for the beaming and ellipsoidal 
variation, integrated over the stellar visible hemisphere 
and wavelength through the Kepler band. In integrating 
the local flux calculated from the EVIL-MC model, we 
assume that the stellar effective temperature is 4550 K, 
and adopt the gravit y darkening exponen t of /3 = 0.093 
for Equation (10) of l.lackso n et al. (1201 2D based on the 
theoretical calculation bv IC1~ i ni 99SD . The free pa¬ 
rameters relevant to the light curve model are a/i?*, the 
transit impact parameter b, limb-darkening coefficients 
ui U 2 and ui — U 2 for the quadratic limb-darkening 
law, q (= Mp/Mi,), Fq, Fi, Acj), the overall normaliza¬ 
tion factor C for the folded light curve, Rp/Ri,, ecosw, 
esinw, and ATc, which represents the small time devia¬ 
tion of the transit center from the ephemeris reported by 
the Kepler team. In our updated ephemeris, the initial 
transit center becomes — ATc. We fix the orbital 
period to be Porb = 6.24668005 d that are derived b y th e 
Kepler team, based on the Q1-Q16 data (see section HTTl) . 

Assuming that the likelihood is proportional to 
exp(—x^/2) in Equation ([1]), we simultaneously model 
the observed radial velocity and light curve, and com¬ 
pute the posterior distribution for each htting parameter 
by implementing Markov Chain Monte Carlo (MCMC) 
simulation. In addition to the above mentioned twelve 
parameters, we add 7 in Equation m to the fitting pa¬ 
rameters. Note that the radial velocity semi-amplitude 
K in Equation ([2]) is related to the planet-to-star mass 


Fp = Fq — Fi 


1 — e sin u) 
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Fig. 3. — Phased radial velocities for Kepler-91 (red points). The 
solid line is the best-fit radial-velocity curve to the data obtained 
by global modeling of photometric and radial-velocity data consid¬ 
ering the velocity trend 7 . The derived 7 (= —0.0612 m s“^day“^) 
was subtracted from the plot. 


ratio q by 


K = 212908.30 


M,/Mq \ 

Porb/day/ 


q sinio 
(1 -b q)i \/l - 


(ms \)t) 


In computing the posteriors, we do not impose priors on 
the fitting parameters except for ui and U 2 ', due to the 
sparse sampling of the long-cadence data 30 minutes) 
and quality of the binned light curve, the limb-darkening 
coefficients are poorly constrained in the absence of pri¬ 
ors, and thus we decide to put Gaussian priors on the 
li mb-darkening coeffic i ents b ased on the theoretical table 
bv iGlaret fc Bloemenl (j2011[ ) as ui -I-M 2 = 0.74±0.01 and 
U 1 — U 2 = 0.49±0.01. In our MCMC algorithm, originally 
developed in iHirano et al.l (I2012D , the step size of each 
fitting parameter is iteratively scaled so that the over¬ 
all acceptance ratio falls between 15% and 35% . After 
running 1 , 000,000 chains, the best-fit value and uncer¬ 
tainty for each fitting parameter are estimated from the 
median, and 15.87 and 84.13 percentiles of the marginal¬ 
ized posterior distribution of that parameter. 

Shortly after we performed the first MGMG trial, we 
noticed a possible trend or drift in the observed radial 
velocities. Figured] plots the best-fit radial velocity curve 
as a function of the orbital phase (p. For clarity, the data 
taken in 2013 are shown in green and those by the 2014 
campaign are plotted in red. While we have detected a 
clear modulation with an amplitude of ^ 50 m s“^, radial 
velocities take lower values for the 2014 data, which is 
evident in the bottom panel indicating the residual of 
the observed radial velocities from the best-fit model. 

In order to confirm the presence of the radial velocity 
drift, we try to fit the data with an additional parameter: 
the radial velocity drift 7 . Adding the drift term 7 x t 
to the right side of Equation ([2]), we performed again 


TABLE 4 

Orbital Parameters eor Kepler-91 


Parameter 

Value (with trend) 

Value (without trend) 

ajR-k 

Q45 

9 oqq+ 0.046 
z.zoo_Q 043 

b 

U.yU‘i±_Q QQgg 

u.yuoz_Q QQg2 

Ul + U2 

n 707+0.010 
^•<'J<_0.0098 

n 7q7n+0-0100 
U.IOIU_0 00gg 

Ui — U2 

0 493+0 010 
U.‘i»0_Q QJQ 

0 493+0 0099 

U.‘i»0_Q QJQ 

q(= Mp/M^) (10-4) 

^•o^_0.i9 

4 72+^-^® 

Fo (10-5) 

o-oo_o 55 

K ^3+0.57 
0.09_q 5g 

Fi (10-5) 

1 41+0.30 
-^•^-^-0.29 

-^•^^-0.28 

A<f> 

0 394+0 029 

u.o»‘i_Q Q33 

0 395+0 028 
u.oyo_0 032 

Rp/Ri, 

0.022861^0030 

u.u^zy‘i_Q 00030 

e cos u) 

0.02801°;°“58 

n no77+0.0040 
u.uz/ <_0.0039 

esincj 

-0.043l“:°“ 

-0.046l°;°“i 

ATc (10~^ day) 

-1.09l°;|l 

-i-iitg:!! 

7 (m s~^) 

-42.512° 

- 33 . 2 II® 

7 (m s“^ day“^) 

—0 061 9+® *^^®^ 
U.UDiZ_Q QQg2 

- 


296 

396 

BIC 

383 

477 

io (deg) 

67.371“;“ 

67.201“;“ 

e 

U.UDiy_Q QQgg 

0.05351°;°°°? 

OJ (deg) 

-57.2l®-( 

-58.81®;? 

Mp (Mjup) 

0.66±0.06 

0.65±0.06 

Rp (Rjup) 

1.40±0.04 

1.41±0.04 


the MCMC simulation to fit both radial velocity and 
light curve. The derived best-fit parameters and their 
uncertainties are summarized in Table Sj To compare 
between the radial velocity models with and without a 
drift term, we compute Bayesian Information Criteria 
(BIC) for the two cases, which is computed as BIC = 
X^+nln(A‘data), where n is the number of free parameters 
and Adata is the number of data points. In the absence 
of a trend, we obtain BIC = 477, and for the case of the 
radial velocity model with 7 , BIC becomes 383; ABIC 
is much larger than 10 , meaning that the model with 7 
is strongly favored. Therefore, we conclude that a radial 
velocity drift (or trend) is present in our dataset, and 
report the best-fit parameters for the case with 7 as the 
final result. We also tried several periods around Porb 
and found that the Porb — cr gave slightly better results in 
the global fitting than Porb did (Ay^ ~ 7). However, the 
resultant parameters for the two cases are well consistent 
within 0 . 2 cr level. 0 

Figures [2| and [3| plot the folded light curve and radial 
velocities (red points) along with their best-fit models 
(blue line). The constant radial velocity offset and trend 
are both removed from the radial velocity data in Figure 
[3| The bottom panel in each figure shows the residuals 
from the best-fit model. From the posterior distribution 
of the fitting parameters, we also estimate the orbital 
and planetary parameters (e.g., the orbital inclination 
io, planet mass Mp and radius Rp) assuming the stellar 
properties reported by asteroseismology (Tabled]) to be 

io = 67.37^0)65 — d -66 ± 0.06 Mjup, and Rp = 

I.40±0.04 Rjup. These estimates are also shown in Table 

SI 


4.3. KOI-1894 


ATc differs by 1.5cr level. 
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-0.5 0 0.5 1 1.5 

Phase 


Fig. 4.— Phased radial velocities for KOI-1894 (red points). The 
solid line is the best-fit radial-velocity curve to the data obtained by 
global modeling of photometric and radial-velocity data with the 
orbital eccentricity fixed to 0 and the stellar-density prior imposed. 


The radial-velocity variations are less visible for KOI- 
1894 owing to the large radial-velocity error compared to 
the small semi-amplitude. The phase-folded light curve 
also shows a very tiny variation, if any, along the orbital 
phase. To extract possible planetary signals, we simulta¬ 
neously model the radial velocities and light curve as in 
the case of Kepler-91. Since the observed transit depth 
and KOI-1894’s stellar radius suggest that the radius of 
the transiting companion (KOI-1894.01) is no more than 
0.65i?j, the reflected/emitted light from the planet is ex¬ 
pected to_be_yery_small^According to Equations (1) - 
(3) in iShporer et al.l (|201lD , the planetary reflection is 
estimated to be < 1 ppm at the location of KOI-1894.01, 
which is smaller than the expected amplitude of the el¬ 
lipsoidal variation 2 ppm). Visual inspection of the 
binned light curve also suggests the absence of secondary 
eclipse. Thus for simplicity, we here neglect the planet 
flux Fp for KOI-1894, and model the folded light curve 
with the ellipsoidal variation (including Doppler boost¬ 
ing) and transit only. Neglecting Fp also helps to avoid 
being stuck at a local minimum during the optimiza¬ 
tion. 

Fixing the system parameters as = 4992 K, M* = 
1.41M|^, the gravity darkening exponent of /3 = 0.097 
(|ClarefJll998D . and Porb = 5.28 789787 d derived with the 
Q1-Q16 data (see section |4T]), we perform MCMC sim¬ 
ulations for KOI-1894 to infer the posterior distributions 
for the fitting parameters. In the present case, we have 
the eleven free parameters: a/P*, 5, ui + U 2 , ui — U 2 , 
q, C, Rp/Ri,, ecosuj, esinw, ATc, and 7. Again, we as¬ 
sume Gaussian priors for the limb-darkening parameters 
as ui + U 2 = 0.72 ± 0.01, and ui — U 2 = 0.39 ± 0.01 


ffl 

0 ) 

cc 


o 

cc 


1.00005 
1 

0.99995 

0.9999 

0.99985 

0.9998 
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0.9997 

0.99965 

0.9996 

6e-05 

4e-05 

2e-05 

0 

-2e-05 

-4e-05 
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Fig. 5.— Phased light curve for KOI-1894 (red points). The solid 
line is the best-fit light-curve model to the data obtained by global 
modeling of photometric and radial-velocity data with the orbital 
eccentricity fixed to 0 and the stellar-density prior imposed. 


from the table bv iClaret fc BloemenI (|2011[ ). Due to the 
weak radial-velocity signal of KOI-1894.01, we found that 
the global fit to the radial velocities and light curve ex¬ 
hibits a degeneracy between the fitting parameters (e.g., 
a/p* and q). Thus, for KOI-1894, we also impose an 
additional prior on the host star’s density from Table [I] 
(p* = ( 0.0259 ± 0.0055 )pm). assurning a Gaussian distri¬ 
bution (iSeaeer &: Mall6n-Ornelasil2003D . Otherwise, the 
fitting algorithm is exactly the same as for Kepler-91. 
The result of the fit is summarized in TablelS The best- 
fit value for the mass ratio is g = (1.02lg 4 g) x 10“"^, 
implying ^ 2cr detection of the planet. The mass and 
radius ratios are translated as Mp = O.lbtg os ^^Jup and 
Rp = O.Gd^o^^ R-jup assuming the stellar mass and ra¬ 
dius in Tablets Table [S] also shows the result of our fit 
in the absence of the Gaussian prior on p*. As expected, 
the mass constraint is slightly weaker for this case. 

The orbital eccentricity is poorly constrained due to 
the small planetary signal and small number of radial 
velocity data, but we can rule out a high eccentricity 
(e > 0.25) with our current datasets. We also try to 
fit the data with a circular orbit (e = 0), and compare 
between the eccentric and circular cases. As shown in 
Table 0 the circular orbital fit yields a slightly better 
constraint on the mass ratio; q = (1.19to)4i) ^ 10“^ 
with the stellar-density prior, leading to a ~ 2.9 (t de¬ 
tection. The mass and radius ratios are translated as 
Mp = 0.18±0.07 Mjup and Rp = 0.66±0.04 Rjup assum¬ 
ing the stellar mass and radius. Based on the values 
for the best-fit parameters, we compare BIG values for 
the two cases; BIG values of 295 and 285 are obtained for 
the eccentric and circular cases (with the stellar-density 
prior), respectively, so that fitting with a circular orbit is 
favored. The observed radial velocities and whole light 
curve are plotted in Figures |4] and 0 with their best-fit 
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Fig. 6 .— Phased light curve for KOI-1894, zoomed in to show 
transit. The points and line are the same as those in figure [5] 


models for e = 0. We also show zoomed-in versions of 
the transit and phase-curve variation in Figures [ 6 ] and [3 
respectively. In Figure [3 the binned flux data used for 
the fit (250 bins) are plotted by the red crosses, and black 
points with errorbars indicate the flux data binned into 
30 bins. We also tested several periods around Porb and 
found that about Porb — 4 (t gave the minimum value 
in the global fitting (x^ 200 ), suggesting that the true 

period might exist around it. Nonetheless, the resultant 
parameters are well consistent with those for the case of 
Porb within 0 . 2(7 level0 

Our global analysis indicates the detection of KOI- 
1894.01 is still marginal, with only ^ 2 — 3a level, and 
the planet mass seems to be mainly constrained by the 
radial velocity data. In order to see if an independent 
estimate for the planet mass from the light curve alone 
gives a comparable result, we fit the folded light curve 
and estimate system parameters without radial velocity 
data. As a result, we obtain q = (1.27lQ g 2 ) ^ 10“"*, 
consistent with the above result by the global fit. 


5. DISCUSSION AND SUMMARY 

We reported the results of high-precision radial- 
velocity measurements with Subaru/HDS for the two Ke¬ 
pler evolved stars Kepler-91 and KOI-1894, which are 
candidates to host transiting planets. 

Based on the simultaneous modeling of the radial- 
velocity data and Kepler light curves, we independently 
confirmed the planetary nature of Kepler-91b. The ra¬ 
dial velocity semi-amplitud e of 55 m s~^ we de¬ 
rived is consistent with what [Barclay et al.l (I2014D found 
{K = 67 ± 11 m s~^) with ^ la, but incompatible 
with the result bv iLillo-Box et aP (j2014bfl who reported 
AT = 93 ± 17 m with > 2a, although the reason for 
this discrepancy is unknown. We newly detected a drift 
of ^ 20 m s“^ between the radial velocities taken at ^ 
1 -yr interval thanks to our better measurement precision 
and longer period of time of observations compared to 
the previous ones. The persistent radial velocity drift is 
suggestive of a possible presence of another companion 

® ATc differs by ~ 3a level. 



phase 

Fig. 7.— Phase-folded light curve for KOI-1894, zoomed in to 
show the flux modulation by ellipsoidal variation. The red crosses 
are binned flux data used for the fit and black points with errorbars 
are data binned into 30 bins. The solid line is the same as that in 
figure [5] 


(planet) outside of Kepler-91b. Due to the lack of data, 
however, we are not able to pin down the period nor 
mass of the companion. Intensive high-precision radial- 
velocity monitoring of the star will uncover the unseen 
companion and provide hints for formation and evolution 
of the planetary system. 

The estimated parameters for Kepler-9 Ib in Table |4] 
are reaso nably in good agreem ent with the previous 
study bv iLillo-Box et all (|2ni4allbl l except for the mass 
ratio q and scaled semi-major axis a/Ri,. This is likely 
because our radial velocity data show a smaller ampli¬ 
tude and planet mass is estimated to be small. Since 
the flux amplitude due to ellipsoida l variation is approx - 
imately proportional to q{a/Ri,)~^ (|ShDorer et al.l[2?iIT[l . 
a small q resulted in the smaller a/i?* in order to explain 
the observed flux amplitude. 

The westward phase-shift {A(f) > 0) of the flux max¬ 
imum in planetary light (from the substellar point) is 
indicative of the inhomogeneous cloud coverage on the 
plan etary surface. Thi s is consistent with the trend 
that lEsteves et al.l (j20I4fl found, claiming that westward 
phase-shifts are preferentially seen for relatively cool 
close-in planets with equ i librium temperatures of Teq ^ 
2500 K: iLillo-Box et all (l2014all estimated the equilib¬ 
rium temperature of Kepler-91b to be ^ 1920 — 2460 K, 
depending on the assumed heat redistribution parame¬ 
ter. The magnitude of the phase-shift (> 0.35) is, how¬ 
ever, unusually large in our best-fit model compared with 
the previously reported values for other close-in planets 
(jLillo-Box et al.ir2014all . The reason is unknown, but the 
huge coverage of the planet surface illuminated by the 
host staiEB might be responsible. Note that we also de¬ 
tected a dip in the folded light curve (Fig ure [2|) around 
(j) = 0.7 reported in lLillo-Box et al.l (l20I4all . and this dip 
should have more or less affected the fitting result. 

As for KOTI894, we did not detect any statistically 

About 70% of the planet surface is always illuminated by the 
host star due to its huge radius. 
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TABLE 5 

Orbital Parameters for KOI-1894 


Parameter 



Value 



e = 0 with p* prior 

e free with p* prior 

e = 0 without p* prior 

e free without p* prior 

a/Rt 

3 70+0.24 

3-6817^ 


3.2417? 

b 

U.UO^_Q 073 

0-56l7o 


0.68l?:i?5 

Ul + U2 

n 79 j^+ 0.0099 
u. 1 ^‘±_0.010 

n 79Q+0.011 
u.izo_0 009g 

0 724+^-^^^ 
u.iz^_0 010 

n yoyi +0.0079 
v.i^i J-_o.o096 

Ui — U2 

n +0.0098 

U.Oyi_Q QJQ 

0.3908t°°K®^ 

u.oy^_0 0099 

0 391+0-010 
u.oyj._o 0096 

q{= Mp/M^) (10-4) 

1 10+0.43 

1 02+^-^^ 
^•^4^-0.49 

J-.^y_0.51 

0 70+0.47 

G. /U_o.4o 

Fo (10-5) 

0 (fixed) 

0 (fixed) 

0 (fixed) 

0 (fixed) 

Fi (10-5) 

0 (fixed) 

0 (fixed) 

0 (fixed) 

0 (fixed) 

A0 

0 (fixed) 

0 (fixed) 

0 (fixed) 

0 (fixed) 

Rp/R* 

0 Qi 777+0.00041 

u.ui 1 1 t _o.00035 

u.uj.roy_Q QQQ35 

u.ur 1 ± i _o.00042 

0.01790l7?“g 

e cos Lo 

0 (fixed) 

U.UOZ_Q Qyy 

0 (fixed) 

0.04717) 

esinoj 

0 (fixed) 

ombtrolt 

0 (fixed) 

O.lOtols 

ATc (10-5 day) 


-0.53l°;5i 

-0.24175 

-0.571?;? 

7 (m s~^) 

-10.671 

-lO.OltJ 

- 11 . 37 ? 

-6.817 


234 

233 

235 

233 

BIG 

285 

295 

286 

294 

io (deg) 

80477” 

80277 

83.57) 

76 2^^-^ 

e 

0 (fixed) 

G.i^y_o 077 

0 (fixed) 

O.ISOIq 097 

U} (deg) 

0 (fixed) 

61+45 

0 (fixed) 

6c;+’^4 

OQ_70 

Afp (Mjup) 

0.18±0.07 

Q 1 r + 0.07 

^•J-o_o.o8 

0.19±0.08 

O-J-O_0.06 

Rp (Rjup) 

0.66±0.04 

O'.0^_0,o3 

0.63±0.04 

o-OO_0.04 



Stellar Radius [Rsun] 

Fig. 8. — Planetary mass plotted against host star’s radius. Red 
circles and blue squares represent true mass and projected mass, 
respectively. Planetary and stellar data were downloaded from 
NASA Exoplanet Archive, and only systems whose stellar mass 
and radius are both reported are plotted. Median and 3cr upper 
limit value are shown for mass of KOI-1894.01. 


significant radial-velocity variations with our measure¬ 
ment precision of 9-15 m s“^ and the number of data 
points. We excluded the possibility of a grazing tran¬ 
sit by a binary companion for the star and set an upper 
limit on the mass of KOI-1894.01 to be 0.39 Mjup by our 
non-detection of radial-velocity variations with 3cr level. 
However, we detected possible radial-velocity variations 


with a semi-amplitude of ^15 m s“^ in phase with ellip¬ 
soidal variations of the star with 2-3(7 level. Although we 
can not say the detection is statistically significant at this 
stage, it suggests that the KOI-1894.01 could be a sub- 
saturn-mass planet. Figure [5] shows distribution of mass 
of exoplanets currently known plotted against their host 
star’s radius. As seen in the figure, KOI-1894.01 could 
be one of the lowest mass planets ever discovered around 
evolved stars together with Kepler-56b, a super-neptune- 
mass planet (Mp = 0.07 Mjup) detected via TTV (transit 
timing variation) method (iHuber et al.ll2013blirq Actu¬ 
ally the stellar par ameters for KOI-18 91 are similar to 
those of Kepler-56 (jHuber et al.ll2013aL M* = 1.32 Mq, 
A* = 4.23 Rq). Confirmation of KOI-1894.01 is highly 
encouraged in order to uncover such a new population of 
sub-saturn and super-neptune planets around relatively 
massive evolved stars, which have rarely been found so 
far either by radial-velocity surveys or transit ones. 

This research is based on data collected at Subaru Tele¬ 
scope, which is operated by National Astronomical Ob¬ 
servatory of Japan (NAOJ). We are grateful to all the 
staff members of Subaru for their support during the 
observations. This research has made use of the SIM- 
BAD database, operated at CDS, Strasbourg, France. 
T.H. is supported by Japan Society for Promotion of Sci¬ 
ence (JSPS) Fellowship for Research (PD:25-3183). K.M. 
is supported by JSPS Research Fellowships for Young 
Scientists (No. 26-7182) and by the Leading Gradu¬ 
ate Course for Frontiers of Mathematical Sciences and 
Physics. 
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